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Abstract 

Propagation of waves through Cantor-set media is investigated by renormalization-group anal- 
ysis. For specific values of wave numbers, transmission coefficients are shown to be governed by 
the logistic map, and in the chaotic region, they show sensitive dependence on small changes of 
parameters of the system such as the index of refraction. For other values of wave numbers, our 
numerical results suggest that light transmits completely or reflects completely by the Cantor-set 
media Coo. It is also shown that transmission coefficients exhibit a local scaling behavior near 
complete transmission if the complete transmission is achieved at a wave number k = k* with a 
rational k* /tt. The scaling function is obtained analytically by using the Euler's totient function, 
and the local scaling behavior is confirmed numerically. 

PACS numbers: 05.45.-a, 42.25.-p, 61.44.-n 
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I. INTRODUCTION 



Wave properties in 
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fractal e, a s, a a a a a b, llOl . Ill| and quasi-periodic 

171 ] structures in one dimension have been of theoretical and practi- 



cal interest over the past two decades. They are typical examples of self-similar structures, 
and physical properties peculiar to them have been explored. In order to observe effects 
of quasiperiodicity, optical experiments using dielectric multilayers of SiC>2 and TiC>2 films 



were performed for the Fibonacci multilayer 



3,0. 



In these experiments, scaling be- 



haviors of the transmission coefficients were observed, which had been predicted by the 



renormalization-group theory [15| . For fractal structures, optical wave propagation on Can- 
tor multilayers has been studied by several authors 2, 3, 4, 5, 6, 7, {J sj]. Self-similar struc- 
tures of the transmission (or reflection) coefficients were obtained numerically 
Moreover, resonant states of light were studied 8|, |9) , which was motivated by an experiment 
using a three-dimensional fractal cavity called the Menger sponge 181 ] . 

In this paper, we reexamine the optical wave propagation through Cantor sequences on 
the basis of the renormalization-group theory. In addition to the self-similar structures 
mentioned above, we find interesting chaotic behaviors of the transmission coefficients an- 
alytically. It is well known that non-linear dynamical systems showing chaotic behaviors 



often have strange attractors with fractal structures [19|, but we show that the reverse is also 
possible. Namely, chaotic dynamics is obtained from a fractal structure. In the following, 
we find that, for specific values of wave numbers, the transmission coefficients show chaotic 
behaviors governed by the logistic map. For these wave numbers, the transmission coeffi- 
cients are very sensitive to small changes of parameters of the system such as the index of 
refraction. This exotic behavior leads to rapid oscillations of the transmission coefficients as 
functions of the index of refraction, which could be observed in an optical experiment. For 
other wave numbers, our numerical study suggests that light eventually transmits completely 
or reflects completely for the infinite generation of the Cantor sequences Coo. 

We will also find intriguing local scaling behaviors of the transmission coefficients near 
complete transmission, which are distinct from the self-similar structures found in Refs. 
4J . The complete transmission can be regarded as a fixed point of the renormalization-group 
equation, and the local scaling behaviors are obtained on the basis of the renormalization- 
group theory. It will be shown that if complete transmission is achieved at a wave number 
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FIG. 1: The first four generations of the Cantor sequences Cj. The Cantor-set media C^ are 
obtained in the j — * oo limit. 



k = k* with a rational k*/tt, the transmission coefficient around the wave number exhibits 
a local scaling behavior. The scaling relation is determined by the rational number k*/tt 
through the Euler's totient function, and an analytic expression of the scaling function will 
be presented. We will also compare the analytical results with numerical data. 

This paper is organized as follows. In SecJTTl we briefly explain the Cantor sequences 
and the Cantor-set media. We then formulate our problem in terms of renormalizat ion- 
group transformation. In Sec lIII[ we classify wave propagation and find chaotic behaviors 
of transmission coefficients governed by the logistic map for specific wave numbers. Optical 
experiments to observe the chaotic behaviors are also proposed. For other wave numbers, 
complete transmission and complete reflection of light in the Cantor-set media Coo are 
numerically suggested. In SecJTVj scaling behaviors of the transmission coefficients near 
complete transmission are analyzed on the basis of the renormalization-group equation. 
Finally, in SecJV] we summarize our results and discuss a generalization of our analysis for 
generalized Cantor-set media. 



II. TRANSFER-MATRIX METHOD AND RENORMALIZATION-GROUP 
TRANSFORMATION 

In this section, we make a formulation of our problem. The Cantor sequences and the 
Cantor-set media are constructed, and the transfer matrix method is introduced to study 
wave propagation through them. Then the renormalization-group equation is defined in 
terms of the transfer matrix method. 

Let us construct the Cantor-set media first. See FigJTJ The procedure of constructing 
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the Cantor set begins with a line segment with unit length (Co in FigJT]). We regard this as 
substrate A. To obtain the first generation Ci, the line segment is divided into three parts. 
The left and the right segments are substrate A, each of which has length 1/3, and the 
middle part, which has length 1/3, is removed. We regard the removed part as substrate 
B. Then the procedure is repeated for each of remaining line segments A to obtain new 
generations. We call the jth generation of the Cantor sequence as Cj. For Cj we have a 
set of 2 J line segments of substrate A, each of which has length 1/3 3 '. By repeating this 
procedure infinite times, we finally obtain the Cantor-set media, Coo, which are self-similar 
and have fractal dimension log 2/ log 3. In the following, we denote the indices of refraction 
of A and B as n A and tib, respectively, and take n-Q = 1 without loss of generality. 

Now consider wave propagation through the Cantor sequence Cj illustrated in FigfSJ For 
simplicity, we suppose that the incident light is linearly polarized. Here E^\ and 
Ejp denote the incident light, the reflected light and the transmitted light, respectively. No 

(2) (2) 

incoming wave from the right E^ exists, E^ = 0. In order to understand light propagation 
through Cj, let us first consider interfaces of two layers in Fig El The electric field for light 
in layer A is given by 

E = E^ exp[ik A x — ioot] + E A 2 ^ exp[— ik A x — iut], (1) 

where k A = n A k is the wave number of light in substrate A and u is frequency of light, (k is 
a wave number of light in the vacuum.) The electric field in a layer B is given by the same 
expression with the subscript A replaced by B. The boundary condition on the interface at 
the position x — I is given by 

E W e ik A i + Efe- ikAl = E^e ikBl + E^e~ iksl , 

n A (E^ ] e ikAl - E { l ] e~ ikAl ) = E^e^ 1 - £ B 2) e" ifcB '. (2) 
By introducing the following variables 

E + = E^ + E&\ E_ = (£« - E^)/i, (3) 

(T5]) is recast into 

T-\n A kl)T AB T(kl) I E+ j . (4) 




A 
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FIG. 2: Electromagnetic wave propagation through the Cantor sequence C,- (j = 2). 



Here T(8) and 7ab are transfer matrices given by 

/ 



T{8) 



cos 5 — sin 5 
sin 5 cos 5 




7-ab = " , (5) 




which represent light propagation within a layer and across an interface A <— B, respectively. 
In a similar manner, from the boundary condition on the interface at x — I + d, we have 

T-\k(l + d))T BA T(n A k(l + d)) ( E+ | , (6) 

B' \ /A 

where 7ba — ^ab is the transfer matrix representing light propagation across an interface 
B <— A. Combining (j3j) with (jSJ), we obtain 

E+ j =T-\h{l + d))% A T{n A kd)T AB T{kl) I E+ | . (7) 

/ B' V / B 

Here note that, for a layer A with thickness d, the phase 5 is given by 5 = n^kd, and for a 
layer B with thickness d, the phase 6 is given by 5 = kd. 

We are now ready to consider light propagation through Cj in Figf2J Repeating the 
similar procedure above, we obtain the following formula for light propagation through Cj. 

= e- ik Mj(k) I E+ I , (8) 

R V / L 

where K±. is defined by ([3]) for and . The real matrix is obtained 

recursively, 

M j+ i(k) = Mj (k/3) T (jfe/3) .M, (Jfe/3) , (9) 





FIG. 3: Electromagnetic wave propagation across interfaces between two layers A and B. 



with the initial condition M. (k) = 7ba7a(^a^)^ab- Using the initial condition and the 
recursive relation, we can show 

det M^k) = 1, (Mj(k)) u = {M 3 {k)) 22 . (10) 

By eliminating from (jSJ), the ratio of the amplitude of the transmitted light to 
that of the incident light E^ for Cy is obtained as 



E ( x ] 2e~ ik 



-2'ij I Hh, - <■/) 

where a,j(k), bj(k) and Cj(k) are components of M.j{k) 



(11) 



M s {k) = | a f ] bM | , aj {kf - b 3 {k) C] {k) = 1. (12) 
Therefore the transmission coefficient Tj = \E^\ 2 /\E^\ 2 of Cj is given by 

where |M,-(Jfe)| 2 = 2 aj {k) 2 + b^k) 2 + Cj {k) 2 . 

In FigH] we illustrate a typical example of transmission coefficients as functions of the 
wave number for C2, C3, and C4. As is shown clearly, we have a scaling behavior of the 
transmission coefficients: If we multiply the wave number by three as one generation in- 
creases, we have a similar structure in the transmission coefficients. To describe the scaling 



behavior properly, we introduce the rescaled wave number kj = 3% for Cj, then from 
we have 

M j+1 (k j+1 ) = Mj(kj)T(^K)Mj(kj), = 1, (MjikMu = (MM)** ( 14 ) 

which can be regarded as the renormalization- group equation describing the scaling behavior. 
The "renormalized" transmission coefficient Tj(kj) for Cj is given by 

4 

Tj{kj) = \Mjk~W+2- (15) 

III. CHAOTIC PROPAGATION, COMPLETE TRANSMISSION, AND COM- 
PLETE REFLECTION OF LIGHT 

In this section we study the scaling behaviors of the wave propagation by using the 
renormalizat ion-group equation ( fl4l) . Using analytical and numerical methods, we find two 
different scaling behaviors depending on n of the rescaled wave number kj = 3 J «: a) For 
K = (m/2 • 3 q )ir with integers m and q, the renormalization-group equation reduces to the 
logistic map describing a chaotic behavior. The renormalized transmission coefficient Tj(kj) 
is very sensitive to parameters of the system, and it changes drastically as the generation 
increases, b) For the other k, it will be found numerically that the renormalized transmission 
coefficient Tj(kj) eventually flows into either T = 1 or T = as j — > oo. 

A. K = (m/2-3 q )ir 

Let us first suppose that m is even, m = 2n. In this case, T(3- ? k) in (}T4"j) takes a 
unique value T(rwr) for j > q [26], which enables us to solve ([H]) analytically. By defining 
Mj = T(3 j+1 K)Mj, (H3D is recast into 

Afj+i = AT/, (j > q), deWj = 1, Mj = T(rm)M h (16) 

and from the relation \Mj\ 2 = \Aij\ 2 , the renormalized transmission coefficient Tj(kj) is 
rewritten as 



T 2 (k2) 




FIG. 4: Transmission coefficients Tj as functions of the wave number k for C2, C3, and C4 (from 
top to bottom) with ua = 2.0. The range of the wave number k for Cj is < k < 3 J 7r. Note that 
the horizontal axes are rescaled by a factor of three. 

To solve f|T6|) . we rewrite Afj in terms of the Pauli matrices Oi [i = 1, 2, 3), 

J\Tj = Xjl + OLj ■ ° ot-j = (oif' ,iaf' ,af'), (18) 
with real Xj an d of 1 {i = 1, 2, 3). From ffl6l) . we have 

Xj+i = X 2 j + ct 2 j, oc j+1 = 2xjOCj, x) ~ <*) = 1, (j > q), (19) 
where ct 2 = {of'') 2 — {otp) 2 + i a f^) 2 - By eliminating at 2 in ffl9l . the map for Xj is obtained, 

Xj+ i = 2 X 2 j -l, (j>q), (20) 
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FIG. 5: (Color online) Transmission coefficients Tj as functions of ua for n = ir/3. The generations 
of the left panel are j = 2, j = 3, and j = 4. We also show transmission coefficients for j = 5, 
j = 10, and j = 100 in the right panel. We plot T m ; n as a function of ha by the thick red line. 




FIG. 6: (Color online) The value of |xi| as a function of tia for k = ir/3. For ua which gives 
< 1 (solid lines), the transmission coefficients show chaotic behaviors. For ha which gives 
> 1 (dashed lines), the transmission coefficients flow into T = for j — > oo. 

and in terms of the solution Xj of fl20l) , ctj and Tj are represented as 

«, = n(2x*)a,, T i= l-{X* + $xf U ~ q) - (21) 



9 




FIG. 7: (Color online) Transmission coefficients Tj as functions of n\ for k = tt/2. The generations 
of the left panel are j = 2, j = 3, and j = 4. We also show transmission coefficients for j = 5, 
j = 10, and j = 100 in the right panel. We plot T m j n as a function of ua by the thick red line. 



Here A and r] are the constants of motion in (TT9|) . 



' =A, -^ = -^=,7, (j>q). (22) 



(2) (2) ' (2) (2) 

a) a q a] a\ 



For \x q \ < 1, they satisfy A 2 + r? 2 < 1, and for \x q \ > 1, A 2 + r/ 2 > 1. 

Equations (1201) and (121"]) can be solved analytically. For < \x q \ < 1, the solution is 



lo/-./ — i I sin[2j ^cos 1 Xg] 

Xj = cos [2 J q cos Xg , «j = r-f ] ] — <* g , 

sm[cos _i XqJ 



T = l-(A 2 + r? 2 ) m 

3 1- (A 2 + r^ 2 ) cos 2 [2^-9 cos- 1 x,]' 



and for \x q \ > 1, 



-l sinh[2 J q cosh |xJ] 

= cosh[2 J cosh \x q \\, a 3 - = — — — — - — rf-<x q , 

smh|cosh \x q \\ 

Tj = < A2 /" 2 )- 1 , . (24) 

3 (A 2 + 77 2 )cosh 2 [2^cosh Ix^l-l 
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FIG. 8: (Color online) The value of |xo| as a function of tia for k = tt/2. For tia which gives 
|xo| < 1 (solid lines), the transmission coefficients show chaotic behaviors. For tia which gives 
|xo| > 1 (dashed lines), the transmission coefficients flow into T = for j — ► oo. 

For \x q \ < 1, the transmission coefficient Tj shows sensitive dependence on parameters 
of the system. To illustrate this, we show the transmission coefficients for q = 1 and n = 1 
case, (namely k = 7r/3), as functions of tia in Figj5l In this case, Xq=i is given by 

1 [2n A n\ VS ( 1 V. f 2n A 7A 

and the region of ha with < 1 is shown in Figj6l It is found clearly that the transmission 
coefficients are very sensitive to ua in the region where |xi| < 1- 

This sensitivity to parameters of the system can be properly understood by introducing 
a new variable Xj = 1 — x]- in terms of the new variable, ( f20|) reduces to the logistic map 
with r = 4, 

X j+1 = rX j (l-X j ), (26) 



which has been studied extensively in the context of chaos 20|, |2l|, |22j, |23j|. The logistic map 
with r = 4 in the interval < X < 1 is known to be very sensitive to the initial condition 
of the system, which implies that Tj for \x q \ < 1 also has the same chaotic property. 

A similar chaotic behavior of the transmission coefficient also appears for odd m — 2n+l, 

K = 7T. (27) 

The matrix T(3 J V) in ([H} is now (-l) j - q T((n + 1/2)tt) for j > q, so A/}(= T(3 j+1 K)Mj) 
becomes Afj = {—\y +1 ^ q T{{n J r l/2)ir)A4j. For this A/}, we have the same equation as (fTfif) . 

A^- +1 = Nj, detA/} = ( 28 ) 
11 



and the same expression of the transmission coefficient as (ITT]) . 

T > = \WT2- < 29) 

By using \ji a 3-> A an d V defined in (TT8I) and (122]) . the same solutions (123]) and (124]) are 
obtained. Thus we have the same class of chaotic behavior. As an example, we show the 
transmission coefficients Tj for k = n/2 (n = q = in (ETJ) as functions of n\ in FigUJ For 
K = n/2, x q =o is given by 

» = 5 (" A + shl (T) ■ (30) 

and the region with \xq\ < 1 is shown in FigJS] We find again that Tj is very sensitive to 
in the region with \xo\ < 1. 

In the chaotic region, |% g | < 1, the transmission coefficient Tj has the lower bound T min 

T min = l-(A 2 + ?? 2 ). (31) 

For k = 7r/3, A and 77 are given by 

^ _ Q _ -v^C^a ~ n A 2 ) sin 2 (7m A /3) + (n A - n£) sin(27m A /3) 

\/3(nA + ^-a 1 ) 2 sin 2 (7mA/3) — (tia + n^ 1 ) sin(27m A /3) — 2y/3 

and for k = 7r/2, 

A = i(n A -^)tan(^), ^ = 0. (33) 

The resultant T min 's as functions of n A are also depicted in Figs J5] and [TJ, respectively. 

From the analytic solution (liMl) . we notice that outside the chaotic region, \x q \ > 1, Tj 
decreases monotonically as the generation j increases. Thus we have complete reflection 
Tj = for j — ► oo. 



B. Other k 



For k 7^ (m/2 • 3 9 )7r, we numerically find that the transmission coefficient Tj flows into 
either T = 1 or T = as the generation j increases. A typical example of our numerical 
results is shown in Figj9] It is also found that if complete transmission of light occurs at a 
certain generation, light always transmits completely in the following generations. To see 
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FIG. 9: (Color online) Transmission coefficients Tj as functions of for n = 7r/4. The generations 
are j = 5, j = 10, and j = 100 (from top to bottom). 

this, let us use the component representation of JAj(kj) as ( jT21 . 

MM = [ % 6j ) , (34) 

which satisfies — bjCj = 1. If Tj = 1, 2a| + + — 2 from (jl~5l) . thus .Mj can be 
parameterized as 

(cos 5j — sin 5j \ 
, aj = cos5j, 6j = — Cj = — sin6j. (35) 
sin 5j cos 5j I 

Substituting this into f[T4"l) yields 

/ cos(25,' + kj) — sin(2<5,' + fc,-) , 

M j+1 (k j+1 ) =[ ) J J [ y | , (36) 

\ sin(25j + fcj) cos(25j + kj) 

which implies T J+ i = 1. 
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FIG. 10: (Color online) Transmission coefficients Tj as functions of tt-a for k = ir/3 + 10 4 . The 
generations are j = 5, j = 10, and j = 100 (from top to bottom). 

When k is near k = (m/2 ■ 3 q )iT, the transmission coefficient is found to show transition 
from the chaotic behavior presented in the previous subsection: The chaotic behavior appears 
in first few generations, but finally the transmission coefficient flows into either complete 
transmission or complete reflection. For example, the transmission coefficient Tj for k = 
7r/3 + 10" 4 is shown in FigJTUl In the fifth generation (j = 5), the transmission coefficient 
shows a chaotic behavior indistinguishable from that for k = 7r/3 (Figj5]), however, in the 
tenth generation (j = 10), it shows a different behavior from that for k = tt/3, then finally 
in the hundredth generation (j = 100), it flows into either T = 1 or T = 0. 

C. Optical experiments 



The chaotic behaviors of the transmission coefficients presented in Sec JIII Al can be ob- 
served experimentally. As was shown in FigsfS] and [TJ the chaotic characteristics in the 
propagation are already evident in the first few generations of the Cantor multilayers, where 
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the transmission coefficients oscillate rapidly as functions of ua- (This oscillation reflects 
the stretching and the folding process of the logistic map.) Therefore, the finite generations 
of the Cantor multilayers are sufficient to observe the chaotic behaviors. Experimentally, 
the finite generations of Cj can be fabricated by using the vacuum deposition on a glass 



substrate 



16 



17| . In a similar manner as the Fibonacci multilayers [16l . Il7j , the Cantor 



multilayers are prepared from Si02 and Ti02 films. NasAlFg and ZnS are also utilized for 
the fabrication of Cj [5]. Although the wave number k needs to satisfy the condition k = 3 j k 
with k = (m/2 • 3 q )ir (m and q are integers) in order to observe the chaotic behaviors on the 
jth generation of the Cantor multilayers, this condition can be met by using tunable lasers. 

From the argument in Sec JIII A| the indices of the refraction of A and B should be set in 
the region with < \x q \ < 1- Suppose that the index of refraction of B is 1, then \ q is given 
by (see Eq.pfl)) 

Xq = ~Tr(T(3« +1 «;)Ag, (37) 

where M. q {= Ai q (k q )) is the solution of Eq. ffT^l) . For example, in the case of k — n/3 
(k = 7r/2), Xq=i {Xq=o) as a function of the index of the refraction of A, tia, is given by 
Eq. fl25l) (Eq. fl30l) ). As was illustrated in Figj6] (FigJHJ), the condition < \x q \ < 1 can be 
met in a broad region of without fine tuning. Therefore, the detection of the chaotic 
behaviors is feasible for the current optical experiments. 

The most impressive chaotic behavior is obtained if satisfies \\ q \ ~ 1- Near n\ 
satisfying \x q \ = 1, all values of T appear in a very narrow region of ha- In particular, 
complete transmission and (almost) complete reflection are nearby in the narrow region. By 
controlling tia, the Cantor multilayers with \x q \ ~ 1 could be used as fast switching devices. 



IV. LOCAL SCALING BEHAVIOR OF TRANSMISSION COEFFICIENTS NEAR 
COMPLETE TRANSMISSION 

In the previous section, we found two classes of behaviors of the transmission coefficients 
where they remain finite: chaotic behaviors and complete transmission. In this section, we 
focus on a behavior of the transmission coefficient near complete transmission. As we showed 
in the previous section, if we have complete transmission of light at a certain generation, 
light always transmits completely after the generation. Thus we can consider that complete 
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transmission is a fixed point of the renormalization-group equation fHM . According to the 
renormalization-group theory, the existence of a fixed point implies that scalings are found 
around it 24] . which turns out to be also true in our problem. In the following, on the 
basis of the renormalization-group equation, it will be shown that if complete transmission 
is achieved at a wave number k* with a rational k*/vt, then the transmission coefficient Tj 
around the wave number exhibits a local scaling behavior which is distinct from the global 
scaling illustrated in FigJH Moreover, we will present the analytic expression of the scaling 
function. 

In order to analyze the scaling behavior of Tj near T = 1, it is convenient to introduce 
new variables Xj = aj(kj), yj = (bj(kj) —Cj(kj))/2 and Zj = (bj(kj) + Cj(kj))/2 for the matrix 
elements of Aij(kj), 

Mfo) bM A. kj :, f , (38) 

V c A k o) a A k i) J 

Since aj(kj), bj(kj) and Cj(kj) satisfy aj(kj) 2 — bj(kj)cj(kj) = 1, the new variables (xj,yj,Zj) 
are constrained on the manifold + Vj = 1 + ^1 shown in FigJTTl By rewriting Xj and yj as 



Z? + 1 cos ipj, yj = J z 2 + 1 sin <pj, (39) 



the map (IT4j) is recast into 



Zj + i = 2Zj\ lz 2 A + 1 cos((fj — &K) 



and the transmission coefficient Tj in (1151) is rewritten as 

Tj = ^f+yj = T+zj' (41) 

The complete transmission Tj = 1 is achieved when Zj = 0. From (f4"0l) . it can be found that 
this occurs either if 1) z\ — with ^ for an integer I (1 < I < j), or 2) Zq — 0. 

First consider the case 1) in detail. Suppose that z\ = with 7^ is realized for 
k = k*. Then from the first equation of ( l4*0i) . it is found that the case 1) is possible only if 
cos(</9;_i — 3' _1 k^) = 0. Thus (fi-i is given by 

, t . 2n + 1 . . 

= 3'-^; + — ^— 7T, (42) 
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FIG. 11: (Color online) Orbits of (|14p are constrained on the manifold x 2 + y 2 = z 1 + 1. Complete 
transmission T = 1 corresponds to z = 0. 



with an integer n. Then using f l40l) with j = I — 1 and z; = 0, we have 

= - e i3! ~ x < (43) 

Hence ^ is obtained as 

^ = 3' -1 4 + 7T, (mod 2tt). (44) 

To determine ^ for j > I, we use the second equation of (j4"D]) . Since Zj = for j > I, the 
second equation of (1401) becomes 

which determines (fj for j > I completely as, 

<pj = -(3 j - 2 j - l 3 l )K* + 2 j - l ipt (j > I; mod 2vr), (46) 



where (pi is given by ff44l) . 

In a similar manner, we can also solve tpj in the case 2). Suppose that z = for k = Kq. 
From the following explicit form of A4o, 

I costiak. — (l/n A ) sinnA/t \ 
Mo(fco) = r BA T A (n A k )T AB =\ A W A; , (47) 

\ nAsinriA/t cosn A K I 

we obtain that zo = ( n A — l/ n A) sin(nA«)/2. Therefore, Kq is given by = mir jn A with an 
integer m. (Note that the indices of the refraction of A and B are different from each other, 
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i.e. % 7^ 1.) Substituting this into (|47|) . we find that xq = cos mix and yo = 0. Thus fo is 
given by 

ipo = inn, (mod 27r). (48) 

Since Zj = for j > 0, we have 

from Eq. fj4*0|) . which determines (fj as 

ipj = ~(3 J ' - 2 J >* + 2Vo, (mod 2tt), (50) 

with (fo = ran. 

Now we study behavior of the transmission coefficient Tj near T = 1. For k near k*, z\ 
becomes nonzero but remains small, so we can neglect higher order terms of zj (i > I) in 
the map (f4"0l) . Up to the next leading order, the map ( HUl) is approximated by 



z j+1 = 2 Zj cos(^- - 3>'k), e^ +1 = e i(2 ^- 3 H (51) 
The latter equation in (loTi) can be solved easily, 

<Pi = -(3 j - 2 j ~ l 3 l )K + 2 j - l <p h (mod 2tt), (52) 

where <pi is given by 

j S^ + TT + Oizt), for/^0 
<£z = < . (mod 27r). (53) 

I rmr + O(z ), for / = 

Here we have determined <pi in a similar manner as (jUJ) and (JUJ), but O(zi) corrections 
appear since Z[ ^ near «;*. Substituting this into the first equation in (|5T|) . we obtain 

Zj +1 — 2zj cos 7® (k) (54) 

where %•(/«) denotes 

(2 ■ 3 J ' - 2^3% - 2^(3^; + 7T + 0(*j)), for I ^ 
7 = ^ ' 1 1 1 (mod tt). (55) 

(2-3 J ' -2- ? >-2-''(m7r + 0(2o)), for I = 
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Therefore, for k near k*, Zj is given by 

i-J 

Zj = z t 



n(2cos T f(«)). (56) 



Using this, we obtain the following transmission coefficient Tj near kJ, 

3-1 

T j = l-zfl[(2cos 1 ?\K)y. (57) 

i=i 

As is proved in Appendix [Aj we can show that if k^/vt is a rational number s/t with 
coprime integers s and £, the phase 7P(/c*) becomes periodic for a sufficiently large i, 

72p(«i)=7{°(«i), (^>r;mod7r), (58) 

where the minimal period p is given by a divisor of ip(t). (Here (p(x) is the Euler's totient 
function [25].) Using this, we obtain a local scaling behavior near complete transmission. To 
see this, rewrite Eq. (15711 by using k = k* + 9 {6 <C 1), 

W) = 1 - zffl (2cos[ 7 f ) «) + iP'KW]) 2 - (59) 
i=i 

If k*/7t is a rational number s/t, then from the relation (l58l) . the product in the right-hand 
side of (15"9"1) for j = pn + m — 1 (m = 1, 2, • • • ,p) is rewritten as 

pn+rra— 2 

n 



J] (2cos[ 7 ? ) («?)+7f («?)*] 

pn+m— 2 ^p{n—l)+m—2 

J] (2cos[ 7 f (*D + 7 f K)^]) II (2«*tf(K» +7f («?)*] 

i=p(n— l)+m— 1 i=i 

p— 1 p(n— l)+m— 2 

n ( 2 ^s[7 P ;Li )+m -i + j^) +7S ) ;_i )+m -i +Q (^)^]) n (^os^vn («n*f 

a=0 i=l 

p— 1 2 p(n— l)+m— 2 

II (2cos[ 7 2 0(m)+m _ 1+a («r) + 7^_ 1)+m _ 1+a «)^]) II (zcob^K) +7f 



a=0 i=Z 

where n (m) is the minimal integer satisfying pn + m — 1 > r. Therefore, we have 

p-i 



(60) 



1 - T pn+m _!(^) = JJ ( 2 COS [7So(m)+m-l+a«) + 7p(Ll)+m-l+a«) ]) ( X ~ T P (n-l)+m-l(#)) • 

(61) 



Q = 
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FIG. 12: (Color online) Transmission coefficients Tj as functions of k for tia, = 2.0 around Kq = ir. 
The generations are j = 3, j = 4, and j = 5. 

From Eq. ([551), 7p(Li)+m-i+a( K *) behaves as 

7$;- 1)+m - 1+ >n - 2 ■ ^H- 1 *", (62) 
for n 1. Thus defining p m \9) as 

/<»>(0) = (2 eos[ 7 « 0(m)+m _ 1+ >n + 2 ■ T6])\ (63) 
0=0 

we obtain 

1 " T pn+m _i(9) = f(m) i3 P(n~l) + m-l e) ^ _ T p ( n _ 1 ) +m _ 1 (0)) (n » 1). (64) 

By renormalizing T pn+m _i(8) as 

T m+m „ x {9) = Tpn+^iCe/^" 1 ^- 1 ), (65) 

Eq.( l64j) is rewritten as 

i - f pn+m _i(fl) = / (m) (^)(i - f^v+m-iie/y)). (66) 

This equation clearly indicates that the ratio (1 — T pn+m _i(8)) / (1 — T p ( n _i) +m _i (#/3 p )) does 
not depend on the generation of Cantor sequences and it has a scaling behavior with the 
scaling function f^ m '(9). 

To illustrate the local scaling behavior of Tj near T = 1 obtained above, we compare 
our formula f|66|) with numerical results for various k*. In Figsjl2] and [13l we show the 

20 



-It It-It II -It II 

e e e 

(b) .j=4 j=9 ,j=14 

4 

3 

1 





FIG. 13: (Color online) (a) Transmission coefficients Tj as functions of 8 for ha = 2.0 around Kq = 
tt. The generations are j = 4, j = 9, and j = 14. (b) The ratio Rj(0) = (l-f i (6»))/(l-f i _i(6'/3)) 
for n\ = 2.0 around Kq = tt. The generations are j = 4, j = 9, and j = 14. The scaling function 
f^ 1 ' (6) is also plotted by the line. 

transmission coefficient Tj for ua = 2 and Kq = it. The transmission coefficients Tj as 
functions of k are presented in Fig {12] and the renormalized one Tj(6) is in Figjl~3l In this 
case, the period p in fl58l) is p = y?(l) = 1, and the scaling function p^(0) is 

/W(0) =4cos 2 (2#). (67) 

For j > 9 we find an excellent agreement between our formula (166]) and numerical data in 
FigfH 

In Figs JT4lfT6l we also present local scaling behaviors for = t and Kq = n/t with 
integers t = 2,3,5. For t = 2,3,5, the period p is given by p = 1,2,4, respectively. The 
scaling functions f^ m \0) are given by 

f {l \6) =4cos 2 (2#), (68) 




- Tt Tt 



e 
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FIG. 14: (Color online) (a) Transmission coefficients Tj as functions of 9 for nx = 2.0 around 
Kq = it/2. The generations are j = 9 and j = 10. (b) The ratio Rj(9) = (l-fj(0))/(l-^-_i(0/3)) 
for «a = 2.0 around Kq = tt/2. The generations are j = 9 and j = 10. The scaling function /^(0) 
is also plotted by the line. 



for t = 2, and 

/W(0) = 16 cos 2 (tt/3 + Q6) cos 2 (2tt/3 + 20), 

/( 2 )(0) = 16 cos 2 (2tt/3 + 60) cos 2 (tt/3 + 26), (69) 

for t = 3, and 

/M(0) = 256 cos 2 (tt/5 + 540) cos 2 (4vr/5 + 180) cos 2 (4vr/5 + 60) cos 2 (vr/5 + 20) , 
/( 2 )(0) = 256 cos 2 (tt/5 + 540) cos 2 (vr/5 + 180) cos 2 (4tt/5 + 60) cos 2 (4tt/5 + 20) , 
/( 3 )(0) = 256 cos 2 (4tt/5 + 540) cos 2 (tt/5 + 180) cos 2 (tt/5 + 60) cos 2 (4vr/5 + 20) , 
/( 4 )(0) = 256 cos 2 (4tt/5 + 540) cos 2 (4tt/5 + 180) cos 2 (tt/5 + 60) cos 2 (vr/5 + 20) , (70) 

for t — 5. Again, we have excellent agreements between the numerical data and our analytical 
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FIG. 15: (Color online) (a) Transmission coefficients Tj as functions of 9 for = 3.0 around Kq = 
tt/3. The generations are j = 14 and j = 16. (b) The ratio Rj(6) = (1 - T i (6»))/(1 - f^sl^/S 2 )) 
for «a = 3.0 around Kq = n/3. The generations are j = 14 and j = 16. The scaling function 
f^(0) is also plotted by the line, (c) Tj for the generations j = 13 and j = 15. (d) Rj for the 
generations j = 13 and j = 15. The scaling function is f^ 2 \9) in this case. 
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71/81 



FIG. 16: (Color online) (a) Transmission coefficients Tj as functions of 9 for = 5.0 around Kq = 
vr/5. The generations are j = 13 and j = 17. (b) The ratio Rj(9) = (1 - T i (6»))/(1 - f i _ 4 (6»/3 4 )) 
for «a = 5.0 around Kq = it/ 5. The generations are j = 13 and j = 17. The scaling function 
f( 2 '(9) is also plotted by the line, (c) Tj for the generations j = 14 and j = 18. (d) Rj for the 
generations j = 14 and j = 18. The scaling function is f^(9) in this case. 
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FIG. 17: (Color online) (a) Transmission coefficients Tj as functions of 6 for n\ = 1/2 around 
k\ = vr. The generations are j = 14 and j = 15. (b) The ratio Rj(9) = (1 - fj(0))/(l - :Z)_i(0/3)) 
for «a = 1/2 around k* = 7r. The generations are j = 14 and j = 15. The scaling function f^(0) 
is also plotted by the line. 



results. 

As an example with / ^ 0, we show a local scaling behavior for tia = 1/2 and k* = it in 
FigHTl Here p = 1 and the scaling function is given by 

f {l \6) =4cos 2 (2#). (71) 

It also reproduces the numerical data excellently. 

We also confirm numerically that if complete transmission is achieved for a wave number 
with an irrational k,*/tt, no scaling behavior is obtained. The case of tia = V% is shown in 
Fig US Although complete transmission is realized at k = -^tt, no scaling behavior of li- 
near the complete transmission is found in FigJTHl 
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FIG. 18: Transmission coefficients Tj as functions of 6 for nx = V% around k* q = ^vr. The 
generations are from j = 4 to j = 9. The horizontal axes are rescaled by a factor of three as a 
generation increases. 

V. SUMMARY AND DISCUSSIONS 

We investigate wave propagation through Cantor-set media on the basis of the 
renormalization-group equation. We analytically find that, for specific wave numbers, the 
transmission coefficients are governed by the logistic map. Especially, in the chaotic region, 
the transmission coefficients show sensitive dependence on small changes of parameters of 
the system such as the index of refraction. For wave numbers near the values giving the 
chaotic behavior, the transmission coefficients again show chaotic behaviors in the first few 
generations. For other values of wave numbers, our numerical results suggest that light 
transmits completely or reflects completely by the Cantor-set media. We also show that the 
transmission coefficients exhibit local scaling behaviors near complete transmission if the 
complete transmission is achieved at a wave number k* with a rational k* /it. The analytic 
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form of the scaling function is determined by k* through the Euler's totient function. 

Finally, we would like to mention that our approach developed here can be extended to 
the generalized Cantor sequences, where the length of each A becomes 1/(5 as one increases 
the generation. For a positive integer j3, in a similar manner as Sec JIII A| it is found that 
the initial wave numbers which give chaotic behaviors are 

2m 



K 



-7r, (m, q : integers). (72) 



G0-l)G0-2)/*» 

In addition, for an odd positive integer /3, chaotic behaviors are also found to appear for 

2m +1 . 

K = FiiF^"' ^■■^^)- ( 73 ) 



APPENDIX A: PERIODICITY OF ~/?\k*) 

In this Appendix, we prove that 7$ («*) has the periodicity (|58|) if and only if k*/tx is 
a rational number. From the explicit forms of 7^ given by fl55|) . it is immediately found 
that if 7} (/c*) has the periodicity (1351) . then /c*/7r should be a rational number. Therefore, 
we will show in the following that if k*/7t is a rational number, then the periodicity ( |58l) is 
obtained. 

To prove this, we use the Euler's theorem: 

N^ M) = 1, (modM), (Al) 
where N and M are mutually prime integers, and <p{M) is the Euler's totient function which 
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counts the number of positive integers not greater than and coprime to M 
totient function satisfies 

(p(nm) = (p(n)ip(m), (A2) 

for coprime positive integers n and m. 

Since O(zi) corrections in Eq. fl55l) disappear for k = 7j- (ft*) is given by 

, (2 • 3* - 2 { - l+2 ■ a 1 ' 1 )*,?, for I ^ 
7?«)= . . !l (modvr). (A3) 

(2 ■ 3 l - 2 J )4, for I = 
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Thus for a rational K*/ir = s/t with coprime integers s and t, we have 



7f(*» = < > ' C 7 ' (mod*). (A4) 



(2f - 2 '-' + ;- 3 '- 1 ) S7r, for Z ^ 



S7T, for / = 

Now decompose t into t = 2 v 3^u where rj, £ and u are integers and it is coprime to 2 and 3. 
Then we obtain 

it, » \ [wk ~ SSS?) S7r > for ^ , , , 

7i U («i)=< > 3i _ € 2l _,x 7 (modvr). (A5) 

Since 3 and 2 r, u, and 2 and 3^u are mutually prime integers, respectively, we have from the 
Euler's theorem 

3 v(a"«) = 1; ( mo d 2^), 

2^ (3eu) = 1, (mod3*u). (A6) 
Moreover, using the relation (1A2I) . we find 

3 ^3««) = (3 ^«))^) = 1? ( mod2 ^), 

2 ^3f«) = (2 ^ u)) ^) = 1? ( mod3 ^ ); ( A7 ) 

namely 

= 1 + 2"uM, 2 v(t) = 1 + 3 5 wA, (A8) 
with integers M and N. From flA8j) . we have 

7 i+v ( t) («r) = ^^FT^ 3S -/+i u J S7r 

= 1? K) + (2 • 3*- ? M - 2 i - /+2 - r '3 / - 1 A)s7r, (A9) 

for / 7^ 0, and 



(0 



f 3 i_ £ 2 i ~ T, \ 
= —— - — - )s7c + (2 • 3 4 ^M - 2 l ""iV) S 7r 
V2'?- 1 m 3% J v y 

= 7i°(«D + ( 2 • 3 " 5m - ^' v N)stt, (A10) 
for / = 0. Therefore, we obtain Eq.( l58l) with p = (p(t) for a sufficiently large i. Here note 
that <f(t) is not the minimal period of 7} (/«*) in general. Thus the minimal period p is a 
divisor of 
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